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The breaking of chiral symmetry of light quarks at zero temperature in presence of strong 
quantizing magnetic field is studied using Nambu-Jona-Lasinio (NJL) model with Thomas-Fermi 
type semi- classical formalism. It is found that the dynamically generated light quark mass can 
never become zero if the Landau levels are populated and increases with the increase of magnetic 
field strength. 

The theoretical investigation of properties of compact stellar objects in presence of strong quantizing magnetic field 
(«~^ ' have gotten a new life after the recent discovery of a few magnetars [1-4]. These stellar objects are believed to be 

■ strongly magnetized young neutron stars. The surface magnetic fields are observed to be > lO^^G. Then it is quite 
P\| : possible that the fields at the core region may go up to lO^^G. The exact source of this strong magnetic field is of 

course yet to be known. These objects are also supposed to be the possible sources of anomalous X-ray and soft 
, gamma emissions (AXP and SGR). If the magnetic field is really so strong, in particular at the core region, they 
^/^ ' must affect most of the important physical properties of such stellar objects and also the rates / cross-sections of 
, elementary processes, e.g., weak and electromagnetic decays / reactions taking place at the core region should also 
I ■ change significantly. 

[ The strong magnetic field affects the equation of state of dense neutron star matter. As a consequence the gross- 
es) ■ properties of neutron stars [5-8], e.g., mass-radius relation, moment of inertia, rotational frequency etc. should change 
^ ' significantly. In the case of compact neutron stars, the phase transition from neutron matter to quark matter which 
may occur at the core region is also affected by strong quantizing magnetic field. It has been shown that a first 
order phase transition initiated by the nucleation of quark matter droplets is absolutely forbidden if the magnetic 
field strength - lO^^G at the core region [9,10]. However, a second order phase transition is allowed, provided the 

■ magnetic field strength < 10^°G. This is of course too high to achieve at the core region. 

I The elementary processes, in particular, the weak and the electromagnetic decays/reactions taking place at the 
> core region of a neutron star are strongly affected by such ultra-strong magnetic fields [11,12]. Since the cooling of 
neutron stars are mainly controlled by neutrino/anti-neutrino emission, the presence of strong quantizing magnetic 
field should affect the thermal history of strongly magnetized neutron stars. Further, the electrical conductivity of 
I neutron star matter which directly controls the evolution of neutron star magnetic field will also change significantly 

' Similar to the study of quark-hadron deconfinement transition inside neutron star core in presence of strong quan- 
] tizing magnetic field, some investigations have also been done on the effect of ultra-strong magnetic field on chiral 
symmetry breaking. In those studies, quantum field theoretic formalism were mainly used [13-19]. In particula, in 
the reference [20], chiral symmetry violation is studied with NJL model with quantum field theoretic approach in 
presence of strong magnetic field and gravitational field and the calculations were done exactly. 
$^ In this article we shall study the effect of strong quantizing magnetic field on the chiral properties of light quark 

matter system with the help of NJL model following semi-classical Thomas-Fermi type mean field approach. Now in 
NJL model, there is no in-built mechanism of color confinement, however, it can produce two chirally distinct phases- 
appropriate for confined quark matter within the bag (not necessarily tiny hadronic bag) and the matter outside the 
bag. These phases are also known as the Wigner phase and spontaneously broken chiral phase respectively. Therefore, 
if one re-formulates the NJL model in presence of strong quantizing magnetic field, it is quite possible to obtain the 
effect of quantizing magnetic field on these two chirally distinct phases and hence obtain the effect of magnetic field 
on chiral symmetry breaking. Further, it is also possible to obtain bag pressure from the difference of vacuum energy 
densities of these two phases and hence its variation with strong magnetic field. Assuming that the confinement and 
spontaneously broken chiral symmetry are synonymous, Bhaduri et. al. obtained some estimate of bag constant from 
the difference of energy densities [21]. In the present article we shall modify these original calculations of Bhaduri 
et. al. [21] and Providencia et. al. [22] to study the breaking of chiral symmetry of light quarks in presence of strong 
magnetic fields and show that the chiral symmetry always remains broken in presence of strong quantizing magnetic 
field if the Landau levels for quarks are populated. Our study is basically an application of the formalism recently 
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developed to study the equation of state of dense fermionic matter of astrophysical interest in presence of strong 
quantizing magnetic field [23]. 

We start with the density matrix p{x,x'), given by 



p{x,x')= ^ ^{x)^\x')0{A-\p. I) 



(1) 



where ip and tjj^ are respectively the negative energy Dirac spinor and the corresponding adjoint, satisiy the equation 

htp = E_il> (2) 
(and similarly for -0^) with the single particle Hamiltonian 

h = -f5T,.{p-qfA)+f3m (3) 

with 



E = 



a 
a 



(4) 



75 and /3 are the usual Dirac matrices, A is the ultra-violet cut off in the momentum integral over and A is the 
electromagnetic field three vector corresponding to the external constant magnetic field of strength along z-axis. 
Here the light quark mass m is assumed to be generated dynamically. Now in presence of strong quantizing magnetic 
field along ^;-direction, the up and down spin negative energy spinors are therefore given by 



ip{x) 



1 



(T,i) 



where 



,(T) 



exp[i{Ej - pyy - Pzz)]v. 



Pzlu \ 

(£_ - m)h 




(5) 



(6) 



and 



[2E_{E_-m)YI^ 





V {E_-m)I,-i J 



(7) 



where E^ = —{pi + ra^ + 2vqfBmY^'^ = —E^, is the single particle energy eigen value, 1^ = 0, 1,2..., are the Landau 
quantum numbers, qf is the magnitude of the charge carried by /th flavor and 



fBr, 



1/4 



(l/!)l/2 



exp 



B„ 



Py 



(8) 



with is the well known Hermite polynomial of order v, and Ly, Lz are respectively length scales along Y and Z 
directions. Now it can very easily be shown that i/ = state is singly degenerate, whereas all other states are doubly 
degenerate. We now express the density matrix, as the modified version of Wigner transform in presence of strong 
quantizing magnetic field, in the following form: 



p{x,x') = ^p{x,x',Py,Pz,iy)exp[i{{t-t')E_ - {y-y')Py- {z - z')pz}] 



(9) 



where the sum is over the momentum components Py, Pz and the Landau quantum number u. Since the momentum 
variables are continuous, the sum over momentum components will be replaced by the corresponding integrals. Then 
using the negative energy up and down spin Dirac spinors, we have 



p{x,x' ,py,pz,v) 



2E. 



-[E_A- PzlzloA + m-^oA - p_L7y7oS]^(A- | pz 



(10) 



2 



where the matrices A and B are given by 



A = 



( lull 




\ 




(11) 



B 



lu-ill 0\ 



h-Jl 





(12) 



where the primes indicate the functions of x' . Now in the evahiation of vacuum energy, we have noticed that it would 
be more convenient to define a quantity jif, similar to the chemical potential for the /th flavor in a multi-quark 
statistical system in presence of strong quantizing magnetic field. Then it is very easy to write 



A = (W^ 



m -2vqfBm} 



(13) 



Since A > 0, it is also possible to express the upper limit of v, which is the maximum value of Landau quantum 
number of the levels occupied by /th flavor, in terms of /i/, m, q/ and B^ and is given by 



max 



2 2 

— 

2qfBm 



(14) 



where [ ] indicates the nearest integer but less than the actual number. Now to obtain the energy density of the 
vacuum, we consider the NJL (chiral) Hamiltonian, given by 



= ^ js(iMi).{pi - qfA) - -Y.5{x, - x^mmj) - /3(i)75(i)/3(i)75(i)] 



(15) 
(16) 



Assuming the magnetic field Bm along z-direction and is constant, we can choose the gauge A^^ = {0,0,xB,0). The 
energy of the vacuum is then given by 



e« = X] / c^a;itri[{75S.(pi - qfA)}pp^]+ei^^ 



(17) 



where Pp^ is given by eqn.(lO) and ei^^ indicates the interaction term, including the exchange interaction. To evaluate 
the vacuum energy, we first calculate the first term of eqn.(17). This quantity is proportional to the trace defined as 
Tr{ph), can easily be evaluated by using p from eqn.(lO) and the single particle Hamiltonian h from eqn.(16). Now 
using the orthonormality relations for the Hermite polynomials at the time of evaluation of integral over dx and also 
using the anti-commutation relations of 7-matrices, we have the first term at zero temperature 



4°^ = 2iv.E ^E(2-M/^P.f 

f=u,d v=0 "'O 



(18) 



where — + 2i^qfB„i, = 3, the number of colors, and E_ = —E^. 

In the evaluation of all the traces in this paper we have used the following important relation: 

Tr(7''7''AiA2..BiB2..) = Tr{A,A2..B,B2..)g'"' , 



(19) 



Tr(7''7V7"^i^2..i?iS2..) = Tr(AiA2..i?iB2..)(5''V^ -rt"" +5'^V"), 



(20) 
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Tr(product of odd 7s with A and/or B) = etc. The other interesting aspects of A and B matrices are: 

i) h^k^^Ti{AiA2) = - A;i,fc2z)Tr(AiA2) 

ii) h^e^^Tr{BiB2) = fcii.fc2±Tr(BiB2) 

iii) fci^fc2''Tr(^iB2) = A:i^fc2MTr(BiA2) = 

iv) pi^/ei''p2.A;2-Tr(^iB2) ¥^0 = {E,,E,,^ - Pizki^)p2±-k2±TT{AiB2) 

These set of relations are very recently obtained by us [23] . Since 7 matrices are traceless and both A and B matrices 
are diagonal with identical blocks, it is very easy to evaluate the above traces of the product of 7-matrices multiplied 
with any number of A and/or B, from any side with any order. 

To evaluate the interaction term, we first consider the direct part which is proportional to Tr{(ipp^)Tr{(3pp^) and 
it is very easy to show that TriP'j^p) — 0. Then using the orthonormality relations for Hermite polynomials and the 
anti-commutation relations for the 7-matrices, we have the direct term 



Fdir = -4W[V(A,m)]2 



(21) 



where 



V(A,m) 



j=u,a 



B„ 



i/=0 



(2 - d,o) 



dpz 



ipl + miy 



(22) 



where rui, = (m^ + 2vqfBjnY/'^ . 

To evaluate the exchange term, we first calculate Tr{{j3pp^){l3pp^)). Now 



2E_ 



- [E-f3A + pzA-jz +mA- p±B'yy] 



(23) 



Then at the time of integration over dxi and dx2 if one uses the orthonormality relations for Hermite polynomials, 
the above trace is given by 



^ _l_ PlzP2z _|_ 
E1E2 



m 



E1E2 



(24) 



where both Ei and E2 are negative. Then in the energy contribution, after integrating over pi^ and P2z, the first 
term gives 



(25) 



Similarly the contribution from second term is given by 



(2-M(A'+™')'^' 



(26) 



and finally, the third term is gives 



m 



f 



A + (A2+m2)i/2 



nil, 



(27) 



To obtain the next term in the exchange part, we evaluate the trace Tr{{(i')5Pp^){(3^^pp^)), which unlike the direct 

case, gives non-zero contribution. Using the anticommutation relations of 7-matrices and as usual with the help of 
ortonormality relations for Hermite polynomials, we have the above trace 



^ ^ PizP2z ^ 



El E2 El E2 



+ m 



1 

Ti 



1 

E^ 



(28) 



The contribution to the interaction energy will again be obtained if we integrate over piz and p2z- Then the first term 
is given by 



4 



27r2 



(29) 



The second term is given by 



27r2 



(30) 



The third term is given by 



27r2 



A + (A2 + mf, 



2U/2 



and finally the fourth and fifth terms, which are identical, given by 



(31) 



N, 



m 



A + (A2 + m2)i/2 



27r2 



(32) 



Then combining all these terms we finally obtain the vacuum energy density. Since the mass m, which is assumed 
to be same for both u and d quarks, is generated dynamically, we obtain this quantity by minimizing the total vacuum 
energy density with respect to m, i.e., by putting dcv/dm = 0. Simplifying this non-linear equation, we finally get 



where 



dcy 
dm 



p = 


27r2 


E 


E(2 


Q = 


27r2 


E«/^™ 


E(2 

I-* 


R = 


Ne 

27r2 


E ^/^"' 


E(2 



-P + 2gQR = 



2r7?/A 



1 



ml (A2 + m2)i/2 
_2 ^ 



2mX 



X - 



ml (A2 + m2)i/2 



(33) 

(34) 
(35) 
(36) 



with 



X = In 



A+(A2+m2)V2 



(37) 



It is therefore obvious from eqn.(33) that the trivial solution to = is not possible in this particular situation. On 
the other hand in a non-magnetic case, eqn.(33) reduces to the well known gap equation, given by 



TO ~ AgVm 



(38) 



where V is the overall contribution of interaction terms. Hence it is obvious that m = 0, the trivial solution exists 
in this non-magnetic or the conventional scenario, investigated by Bhaduri et. al. [21]. The phase with m = is the 
Wigner phase and m 7^ is the so called Goldstone phase, which further gives 



45V =1 



(39) 



which is nothing but the well known gap equation used in BCS theory. 

The non-cxistcnce of trivial solution indicates the spontaneously broken chiral symmetry in presence of strong 
quantizing magnetic field. Therefore as soon as the Landau levels are populated for light quarks in presence of 
external magnetic field, the chiral symmetry gets broken, the quarks become massive and the mass to (assumed to be 
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same for both u and d quarks) is generated dynamically. As mentioned earlier, that such investigations were done with 
quantum field theory. On the other hand the approach followed by us in this article is in some sense semi-classical 
and have not been done before. 

Therefore we may conclude that the Wigner phase does not exist in the case of relativistic Landau dia-magnetic 
system. Further, if the deconfinement transition and restoration of chiral symmetry occur simultaneously, then it 
puts a big question mark whether the idea of bag model is applicable at all in presence of strong quantizing magnetic 
field. 

To illustrate the variation of dynamical quark mass with magnetic field, we consider the relation 

= < # > (40) 

where m,r is the pion mass, mo is the quark current mass and /,r is the pion decay constant. Using the spinor solutions 
given by eqns.(6) and (7) we get 

2 2mom Nc ^ , 

We have now solved the eqns.(33) and (41) numerically to obtain A and m for various values of magnetic field. We 
have considered the following sets of numerical values for the parameters. The current quark mass rrio = lOMcV, 
pion mass = 140MeV, pion decay constant = 93MeV, coupling constant g = lOGeV"-^ and electron mass 
rUe = 0.5MeV. In fig.(l) we shown the variation of dynamically generated quark mass with the strength of magnetic 
field. As it is evident that the dynamical quark mass never goes to zero and diverges beyond ~ lO^^G. 

1 I 1 1 iiii| — I I 1 1 iiii| — I I 1 1 iiii| — I I 1 1 iiii| — I I 1 1 nil 



A + (A^ + mi 



,1/2 



nil, 



(41) 
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I I 1 1 iiii| — I I 1 1 iiii| — I I 1 1 iiii| — I I 1 1 iiii| — I I 1 1 nil 

1 10 10 ' 10 ^ 10 ' 10 * 

FIG. 1. The variation of dynamically generated quark mass with the strength of magnetic field (expressed in terms of 
= 4.4 X 10"G.) 
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